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Abstract 

In energy-constrained wireless sensor networks, energy efficiency is critical for 
prolonging the network lifetime. A family of ant colony algorithms called 
DAACA for data aggregation are proposed in this paper. DAACA consists of 
three phases: initialization, packets transmissions and operations on pheromones. 
In the transmission phase, each node estimates the remaining energy and the 
amount of pheromones of neighbor nodes to compute the probabilities for dy- 
namically selecting the next hop. After certain rounds of transmissions, the 
pheromones adjustments are performed, which take the advantages of both 
global and local merits for evaporating or depositing pheromones. Four dif- 
ferent pheromones adjustment strategies which constitute DAACA family are 
designed to prolong the network lifetime. Experimental results indicate that, 
compared with other data aggregation algorithms, DAACA shows higher su- 
periority on average degree of nodes, energy efficiency, prolonging the network 
lifetime, computation complexity and success ratio of one hop transmission. At 
last, the features of DAACA are analyzed. 

Keywords: wireless sensor networks, data aggregation, ant colony 
optimization, energy efficiency, network lifetime 



1. Introduction 

The main task of wireless sensor network(WSN) is to detect and report 
the events of the physical world. In most cases, nodes are batteries powered Q 
with limited energy resources. Suppose that when a node runs out of power and 
stops working, the original transmission paths will be changed. Nodes nearby 
will suffer from heavier work, because of sharing responsibility of the exhausted 
node which casts heavy burden of them. The energy dissipating rate of these 
nodes will become faster. This process spreads which will cause the packets 
loss or even network congestion. Moreover, the performance of the network 
depends on the persistence of the sensors to a large extent. Hence, the main 
challenge for the energy-constrained network is to design energy-efficient routing 
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protocols 0, H, 0, H| which guarantees the persistence and balances the energy 
consumption of the network. 

In recent years, methods of data aggregation [t| are attracting the attentions 
of the researchers. A large number of works of data aggregation protocols in 
wireless sensor networks for energy preservation have been published. 

To begin with, we will provide a brief background of the data aggregation 
methods of WSN. The key idea is to combine the data coming from different 
sources which eliminates redundancy 10], minimizes the number of transmis- 
sions and thus saves energy. This method shifts the focus from the traditional 
address centric approaches for networking to a data centric approach. 



Some existing methods [llj, |12|, [l3j, LLJ] focus on establishing and maintaining 



topology of the data aggregation tree in which data packets are transmitted and 
aggregated to specific nodes. In [llj], the authors propose a method which gen- 
erates Minimum Spanning Tree (MST) by implementing Prim algorithm from 
a global view. In [12j, each node locally forms MST to establish the aggrega- 
tion topology to eliminate energy costs. In (l3| . the authors integrate LMST 
with RNGpjJ to establish the network topology. In [14j, the authors put for- 
ward an ant colony based data aggregation method which takes advantages 
of pheromones to heuristically set up the network topology. In the aforemen- 
tioned methods, packets are transmitted in pre-allocated paths and aggregated 
in some specific nodes. Nevertheless, such topologies consume much energy in 
maintenance. Nodes in the aggregation tree are actually facing large energy 
consumption. The common problems of the existing methods are: 

1 . The energy consumptions of constructing and maintaining network topolo- 
gies are high. 

2. Issues of balancing energy consumption (l6 . [l7 1 are not considered. 



The ant colony optimization(ACO) 18] is one of the most useful swarm 
intelligence (l9| which has been successfully applied in many optimization prob- 
lems such as TSP[2(|, CVRP[2l[ as well as routing in wireless sensor networks (2^, 
23l 24 1 . In this paper, a measure of data aggregation based on ant colony opti- 



mization is proposed which overcomes shortcomings mentioned in 0, [2^, 1261127 
H] ■ We expand the ant colony algorithms [2{| by proposing a family of energy- 
efficient Data Aggregation Ant Colony Algorithms (DAACA for short), which 
contains Basic- DAACA(i.e. the Basic algorithm of DAACA), ES-DAACA(i.e. 
Elitist Strategy based DAACA), MM-DAACA(i.e. Max-Min based DAACA) 
and ACS-DAACA(i.e. Ant Colony System based DAACA). 

DAACA include three phases: 1) the initialization, 2) packets transmissions 
and 3) operations on pheromones. Initially, all the sensors are deployed ran- 
domly, then they set up their routing tables in a self-organization manner. We 
consider each node as an artificial ant which is raised in ant colony optimization 
for solving the global optimization problem. The packets transmitted between 
nodes are considered as tools for updating or adjusting the pheromones. Each 
node estimates the remaining energy and the amount of pheromones to compute 
the probabilities for dynamically selecting the next hop. After certain rounds 
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of transmissions, the adjustments of pheromones are performed, which origi- 
nally combines the advantages of global and local adjustments for evaporating 
or depositing pheromones. Four different pheromones adjustment strategies are 
designed to achieve the global optimal in prolonging network lifetime. 
The contributions of the paper can be summarized as follows: 

1. We proposed a family of ant colony algorithms for data aggregation which 
aims at saving energy and prolonging network lifetime. 

2. We devise the evaporating and the depositing pheromones approaches 
which take advantages of global and local merits. 

3. We design three heuristic methods(i.e. ES-DAACA,MM-DAACA,ACS- 
DAACA) for improving the performance of the basic algorithm(i.e. Basic- 
DAACA) which specialize the way of evaporating and depositing pheromones. 

4. We develop a platform to compare the characteristics (e.g. average energy 
cost, network lifetime, average degree and so on) of DAACA family with 
other data aggregation methods(e.g. PEDAP, PEDAP-PA, L-PEDAP and 
so on). 

The remainder of the paper proceeds as follows. Section [5] gives an overview 
of the state of the art of data aggregation in wireless sensor networks. Section[3] 
introduces the preliminaries of the paper. Section U presents the system model 
and defines the problem. Section[5]describes DAACA in details. The simulation 
results and performance analysis are presented in section [6l Section [7] concludes 
the paper. 



2. Related Work 



In literatures, many algorithms are propos ed to prolong the lifetime of the 
network based on data aggregation [3^, [31 , 16|, 11, 12 , 13 , 14 1 . Some references 



are concentrating on establishing energy efficient topologies to save energy. 



In [111 ], an algorithm named PEDAP (Power Efficient Data gathering and 
Aggregation Protocol) is proposed, which generates Minimum Spanning Tree 
(MST) by utilizing Prim algorithm. It is a near optimal minimum spanning 
tree based routing scheme. But it cannot guarantee the energy efficiency of 
the network, then an energy-aware version of PEDAP, namely PEDAP-PA is 
proposed to maximize the network lifetime. However, the sink node needs to 
periodically broadcast and calculate the MST of the network, thereby, the work- 
load of the sink node is high. 

In [l2| , a Local Minimum Spanning Tree algorithm called LMST is presented 
to establish the network topology, although it can effectively reduce the average 
degree of nodes, some prominent problems still emerge. Each node needs to 
periodically calculate and update the MST (Minimum Spanning Tree) locally 
which leads to a high computational overhead for each node. Moreover, each 
node needs to communicate with its neighbors to obtain the energy condition 
of nei ghbo rs, which still costs much energy. 

In [l3| , a localized, self organizing, robust and energy efficient data aggrega- 
tion algorithm named L-PEDAP is proposed which combines LMST with RNG 
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15| . Although it is proved to have the capability of prolonging the network life- 



time, its topology construction procedure is nearly identical with that of LMST, 
hence, it cannot be considered as an energy-efficient algorithm. 

To reduce the cost of constructing and maintaining the energy-efficient topol- 
ogy, a heuristic method called ACA (Ant Colony Algorithm) which is based on 
ant colony algorithm is proposed in 14]. In the ant colony optimization, a 
colony of artificial ants is used to construct solutions guided by the pheromones 



trails and heuristic information [32[. This behavior enables ants to find the 



shortest paths between the food source and the nest in a random manner. The 
author initiatively designs the rules for depositing and evaporating pheromones. 
However, the requirements of depositing pheromones can be easily met. When 
depositing pheromones, each node needs to communicate with its neighbors. 
Therefore there are lots of communications within the pheromones adjusting 
period, which leads to much energy cost overhead. Therefore, it cannot be 
considered as an energy-efficient algorithm either. 

Moreover, although some other methods do not depend on energy efficient 
topologies to save energy, some literatures are still valuable to be mentioned. In 



311 ] . the authors describe the design and implementation of a running system 



for energy-efficient surveillance which allows a group of cooperating sensor to 
detect and track the positions of moving vehicles. That method can trade off 
the energy-awareness and surveillance performance by adjusting the sensitiv- 
ity of the system. In (30| . the data aggregation is motivated by optimizations 
of aggregation queries. The authors present variety of techniques to improve 
the reliability and performance of the proposed solutions. In [16|, the authors 
explore two methods to further reduce energy consumption in the context of 
network aggregation in sensor networks. Firstly, a group-aware network config- 
uration method is designed, which groups the sensors into clusters. Secondly, a 
framework to use temporal coherency tolerances in conjunction with aggregation 
to save energy is proposed. 



3. Preliminaries 

In this section, we illustrate related knowledge in this paper. 

3.1. Ant Colony Optimization 

The idea of Ant Colony Optimization is originated by Marco Dorigo[33j], in 
observing the exploitation of seeking shortest path between their nest to the food 
source although limited cognitive abilities are owned by ants. In Ant Colony Op- 
timization, artificial ants are implemented to search for heuristic solutions of the 
problems. Originally, artificial ants are utilized to find the solutions stochas- 
tically based on the probabilities whose calculations are defined according to 
certain circumstances(i.e. according to the characteristic of the problem). After 
a feasible solution is constructed, the evaporation and releasing of pheromones 
are executed which dynamically update the amount of pheromones. Therefore, 
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the solution will gradually approach to the global optimal solution of the prob- 
lem. This procedure repeats until certain number of ants have finished their 
trails. The process of ACO can be summarized as Algorithm [TJ 

Algorithm 1 ACO Pseudocode 
1: Initialization of data structure and parameters 
2: while terminate conditions do 
3: Construct solutions 
4: Update statistics 
5: Update pheromones 
6: end while 



In general, we can summarize the advantages of ant colony optimization. 

1. ACO is a probabilistic mechanism for solving computational problems 
which can be reduced to finding good paths in graphs. 

2. The solution of ACO will gradually approach to the global optimal solution 
in terms of adjusting the quantities of pheromones. 

3. In ACO, each ant is considered to have individually limited cognitive abil- 
ity. 

4. The process of ACO is the cooperation of ants which gradually find optimal 
solutions in graph. 

In wireless sensor networks, a large number of nodes are deployed. To find 
the global optimal prolong lifetime solution can be transformed into forming 
data aggregation tree by global optimization. Hence, firstly, our aim is to find 
a solution which can approach to global optimization. Secondly, to prolong the 
lifetime of the network, we need to balance the energy cost of each node, there- 
fore, we need to establish dynamical topology of the networks. Thirdly, each 
node has limited computation capability, it only hold the local information(i.e. 
the information of neighbors). What's more, the network can be abstracted as 
a graph where the vertexes stand for the nodes and the edges represent the link 
between nodes. Last but not least, researches of the ACO in wireless sensor 



networks [1J, [3J, |35[ lay a good foundation for our work. From the analysis 
above, we conclude that ACO is suitable for data aggregation in wireless sen- 
sor network. Therefore, we employ ant colony optimization in resolving data 
aggregation in wireless sensor network. 

3.2. Ant Colony System 

Ant Colony System(ACS) applies the mechanisms of ACO but some changes 
has been made to overcome the drawbacks of ACO and enhance the performance 
of ACO. The main differences can be summarized as: 

1. It applies a more active behavior, which can make full use of the accumu- 
lated pheromones in the path. 

2. The evaporating and releasing of the pheromones are only carried out in 
the most optimal path so far. 
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3. When an ant passes through a trail, the pheromone of that trail will be 
reduced which aims at enhancing the possibility of finding more optimal 
solutions in other trails. 

Inspired by the features of ACS, we implement the ant colony system into 
our data aggregation method in section f5. 51 to enhance the performance. 

3.3. AC A Algorithm 



In |14[, the authors propose an data aggregation(ACA). First, each compute 
the hop count to the sink node. When the source node wants to send data, it 
will select the next hop node by the probability which mainly depends on hop 
count of the neighbor nodes to the sink. 

r(i,j) x J]{i,jY 



Pk(i,j) 



S uG Ar i r(i, u) x n(i, u)P 



where r(i,j) is the pheromones level from node i to node j, is the inverse 

of the hop count from node j to the sink node adding one, Nj is the number 
of neighbors of the node i, and /3 is a parameter which determines the relative 
influence of heuristic values rj(i,j), > 0. Then the authors design the message 
communication mechanism which increase the probability of intersection of the 
routing path to promote the results of data aggregation. Finally, the evapo- 
ration and releasing of pheronmones are defined which dynamically adjust the 
quantities of pheromones in the paths. A node who has a lower hop count and 
happen to exist in the path are more likely to get pheromone which enhances 
the probability of obtaining messages in later transmissions. Whereas a node 
cannot receive a message for more and certain round, the pheromones of such 
node will be evaporate, which reduce the possibility of receiving of nodes of such 
nodes. 

To aggregate size of packets, in this paper, we develop a family of ant colony 
based data aggregation algorithm DAACA. The similarities and differences be- 
tween DAACA family and ACO can be summarized as follows: Similarities: 

1. Ant Colony Optimization methods are employed which aims at establish- 
ing data aggregation tree. 

2. Methods of releasing and evaporating of pheromones are developed. 

3. The ultimate targets are to save energy. 

Differences: 

1. DAACA focus on balancing the energy cost of the whole network to pro- 
long the network lifetime whereas ACO does not take the lifetime into 
account. 

2. When selecting the next hop, the remaining energy, the distance between 
the neighbor to the sender are considered which aims at fairly selecting 
the next hop based on energy cost in DAACA, whereas ACA only cares 
about the distance and the quantities of the pheromones. 
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Table 1: Energy Dissipated Table 



Symbol 


Meaning 


Energy Dissipated 


Erx-elec 
E Rx — elec 


Transmitter Electronics 
Receiver Electronics 
Transmit Amplifier 


50nJ/bit 
50nJ/bit 
lOOJ/bit/m 2 



3. In ES-DAACA, MM-DAACA and ACS-DAACA, we take advantages of 
accumulated pheromones to increase the probability of selecting the global 
optimal solution. But ACO does not take any actions to obtain global 
optimal solutions. 

4. In MM-DAACA and ACS-DAACA, the range of the quantities of pheromones 
are constrained which aims at avoiding trapping into local optimal solu- 
tions whereas ACO does not limit the ranges of pheromones. 

5. In ACS-DAACA, we reduce the probabilities of the nodes who locates 
in the aggregation tree to diversify the solutions to get more optimal so- 
lutions^. e. the data aggregation tree). On the contrary, ACO can not 
guarantee the diversity of the solutions. 

3.4- Energy Consumption Model 

Many energy models [27j are used for energy consumption in wireless sensor 
networks. In our work, we employ the model in which a packet with the size k 
to be transmitted in the distance of d will consume: 

E Tx (k,d) = E Tx 

— elec X /C + 6 am p x k x d 2 (1) 

for the transmitter and 

E Rx (k) = E Rx - e le, c x k (2) 
for the receiver respectively. The definitions of the symbols are listed in table 

m 

4. System Model 

We consider a scenario where the sensor nodes are homogenous and energy- 
constrained. All the nodes (including the sink node) are stationary and ran- 
domly distributed. Further, replenishing energy via replacing batteries on hun- 
dreds of nodes (in possibly harsh terrains) is infeasible. The basic operation 
in such a system is the systematic gathering of sensed data to be eventually 
transmitted to a base station for processing. The key challenge in such data 
gathering is conserving the sensor energies, so as to maximize their lifetimes. 

4--1- Definitions and Assumptions 

In this section, some definitions and assumptions are defined below. 
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4- 1.1. Definition 

Remaining Energy: The initial energy minus the energy consumed in 
transmission process. 

The remaining energy reflects the usage of energy. Moreover it can indicate 
when the network exhausts. 

Energy Difference: The maximum remaining energy minus the minimum 
remaining energy in the network. 

The energy difference indicates the balanced usage of energy in the network. 
Small energy difference shows the balanced usage of energy of all nodes in the 
network and vice versa. 

4-1.2. Assumption 

Assumption 1: In one aggregation procedure, n packets with the size of k 
will be eventually compressed into one with the size of k. 

If the packet size is bigger than k, it will be sliced into pieces with constant 
size of k. If the size if smaller than k, it will be enlarged into size of k. This 
assumption is identical with that of [321 ] 




Assumption 2: Only the transmission process costs energy. 

Since we mainly focus on the energy cost in data aggregation to prolong the 
lifetime, therefore we omit effects of the energy costs of nodes of their own(e.g. 
calculating, waiting and so on) . We assume that the energy cost can only happen 
when transmitting or receiving packets. 

4-2. Problem Statement 

The network is modeled as a visibility graph G — (V, E), where V is the set of 
sensor nodes. Each node has the maximum transmission range denoted by R, by 
which it can setup its neighbor set and routing table, represents the distance 
between node i and node j which is smaller than R, meanwhile the neighbor set 
(denoted as Nb(i) ) of node i indicates a set of nodes whose distances are less 
than R(i.e. Nb(i) — {j\eij < R}). The set of comprises the E. Each source 
periodically senses and collects the data from the surroundings, then sends the 
data to the next hop until the data reaches the sink node. The goal is to find 
an energy efficient routing path to prolong the network lifetime. The best way 
to solve this problem is to aggregate the multiple data packets and reduce the 
traffic in the network. Thus, multiple transmission paths should be merged into 
smaller ones to relieve the energy consumptions of data transmissions. Multiple 
sources periodically send data to their neighbors until the data packets reach 
the sink node in a hop by hop manner. 

5. Ant Colony Algorithms for Data Aggregation 

Overall we propose the system architecture of DAACA. Then we design a 
basic algorithm, after that three heuristic algorithms are proposed to enhance 
the performance of the basic algorithm. 
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Since each algorithm in DAACA family has the same algorithm structure, 
differences only exist in the process of evaporating and depositing phcromones. 
We firstly illustrate algorithm structure of DAACA family and then introduce 
each algorithm in details. The symbols used in DAACA are illustrated in table 

m 



Table 2: Related Definitions 



Symbol 


Notations 


R 


The maximum transmission range. 


id 


The identity of the node. 


Nb(i) 


The neighbor nodes of node i. 


round 


The transmission times. 


roundToU pdate 


The number of rounds to update the 




pheromones. 


P(i,j) 


The probability for node i to select node 




j as the next hop. 




The pheromone of node i to node j. 


E 'distanced 3) 


The energy distance from node i to 


node j. 


Edistanceii i j) 


The Euclidian distance between node i 




and node j. 


ei(i) 


The current value of node i divided by 








The estimation energy of node i to node 




j divided by E lmt . 


e ij 


The distance between node i and node 


-&cur \Jr) 


j- 

The current energy of node i. 


-^estimate y"> 3) 


The evaluation energy of node i to node 




j- 

The initial energy of each node. 


r{i,j) 


The inverse value of E' distance {i, j) 


Times{i, j) 


The transmission times from node i to 




node j. 


MpCost 


The minimum energy cost in a path. 


VMax 


Upper bound of phcromones. 


f]M in 


Lower bound of pheromones. 


conj(i) 


The conjunction times of node i. 


source 


The source node(s for short) of the net- 




work. 


destination 


The destination node (d for short) of 




the network. 


current 


The current node(c for short) who re- 




ceives the data packet. 


Continued on next page 
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Table 2 continued from previous page 



Symbol 


Notations 


next 
p 

1 mm 


The next nodc(n for short)which cur- 
rent will send the packet to. 
The minimum energy cost path. 



5.1. DAACA Algorithm Structure 

Similarly in wireless sensor network, nodes are considered as the artificial 
ants. The amount of the pheromones of links are recorded in the routing table, 
whenever a node receives a packet, it will figure out the next hop based on the 
information of the pheromones, the distance and the remaining energy according 
to the routing table. Then each node will update the information. After certain 
rounds of transmissions, procedure of evaporating and depositing pheromones 
will be applied which help DAACA to adjust the amount of the pheromones. 
The objective of this process is to guide the aggregation routing path to approach 
to the global optimal routing path which will conserve energy for the network 
to some extent. Next we will illustrate the realizations of DAACA in details. 

Initially every node constructs its own neighbor set and routing table by 
broadcasting its id (i.e. the identity of the node) and location information 
within R. Afterwards, the transmission process starts. Packets are sent from 
the source nodes to the sink node in each round. When a node receives a packet, 
it will evaluate contents of the routing table (e.g. the remaining energies, the 
amount of pheromones) to calculate the transmission probabilities of selecting 
the next hop. Then one of its neighbors will be selected. When the number 
of round equals to the roundTolfpdate, all nodes will perform the evaporating 
and depositing operations. Each node evaporates the pheromones of its neighbor 
and deposits pheromones according to different specific conditions (referring to 
section [531 - 15.51 in the following sections). After that each node will update the 
pheromones of its neighbors. 

Algorithm [2] illustrates the structure of DAACA. Firstly, the initialization 
of the network is carried out. Each node broadcasts "hello" message to its 
neighbors. Then only the nodes who are nearer to the sink node will be recorded 
in the routing table and the contents of the table are initialized. Then the 
transmission begins, the source nodes periodically send the data packets to the 
sink node. When a relay node receives a packet, it will calculate the probability 
of selecting the next hop based on the elements in the routing table (e.g the 
distance, the amount of the pheromones, estimated remaining energy et.al). The 
transmission proceeds until the data are forwarded to the sink node (i.e. line|3} 
linc !18|) . After transmitting roundToUpdate rounds, adjustments of pheromones 
are carried out (i.e. line SOI line Each node updates the routing table to 
keep the information fresh. At last each node broadcasts its actual energy, by 
which the neighbor nodes will update of the contents of the routing table (i.e. 
lineElhneEl). 
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Algorithm 2 DAACA Structure 
1: Network Initialization 
Node Initialization 
Neighbor Initialization 
Routing Table Initialization 
2: The source node begin to send its data packets to the destination hop-by- 
hop. 

3: for all t e Nbr(s) do 

4: s evaluates the energy of t. 

5: s calculates p(s, t). 

6: s selects the next hop node n based on p(s, t) 
1: end for 

8: s — > nils sends the packets to n 
9: while s ^ d do 
10: c = n 

li: c aggregates data packets 
12: for all t G Nbr(c) do 
13: c evaluates the energy of t. 
14: c calculates p(c, t). 

15: c selects the next hop node n based on p(c, t) 

16: end for 

17: c — > n 

18: end while 

19: round = round + 1 

20: if round = roundTol/pdate then 

21: Evaporating Pheromones. 

22: Depositing Pheromones. 

23: Updating Routing Table. 

24: Energy Broadcasting. 

25: end if 
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-DAACA 



ES-DAACA 





Basic-DAACA 


\ 


ACS-DAACA 











ES-DAACA = Basic-DAACA + global pheromones depositions 
MM=DAACA = ES-DAACA + limitations of pheromones 
ACS-DAACA = Basic-DAACA + global & local pheromones adjustments 



Figure 1: The relationship between Basic-DAACA, ES-DAACA, MM-DAACA, ACS-DAACA 



In the following sections, we present four ant colony heuristic algorithms 
(i.e. Basic-DAACA, ES-DAACA, MM-DAACA, ACS-DAACA) for adjusting 
the pheromones especially the evaporating and depositing of pheromones. Basic- 
DAACA is the basic algorithm, other algorithms are designed based on Basic- 
DAACA. The relationship between Basic-DAACA, ES-DAACA, MM-DAACA 
and ACS-DAACA can be formalized in figure [1] 

ES-DAACA adds additional global pheromones depositions into Basic-DAACA. 
MM-DAACA limits the ranges of pheromones in ES-DAACA. ACS-DAACA 
adds both global and location pheromones adjustments in Basic-DAACA. Fol- 
lowings are explicit illustrations of DAACA family. 



5.2. Basic-DAACA 

Basic-DAACA is the basic algorithm of DAACA family. For the pheromone 
is the most critical part in adjusting the probabilities in the routing table, we 
firstly illustrate the constitution of the routing table. The routing table contains 
the following elements: { id, E' distance (i, j), E estlmate (i, j), T)(i,j), p(i,j) }, j S 
Nb(i). id refers to the identity of the node. p(i,j) is the probability for node i 
to select node j as the next hop which is computed as follows: 

P(M)= iM^hit (3 ) 

S \r{i,j) a x rf[i,j) p ] 

j£Nbr(t) 

a and j3 are two parameters which determine the relative influence of r(i,j) 
and eta(i,j). r(i,j) represents the inverse value of the energy distance E' distance (i, j) 
between node i and j, which holds: 

distance \ ' J t 

E'distance&j) can be computed as: 

ELtanceihj) = , ~ ). H O < Bl < 1, < e 2 < 1 5 

Where: 

/ .n E cur (i) 
ei(*) = -=; (6) 
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e*(t,j) = = (7) 

^init 

Edistance(i, j) is the energy distance which can be calculated as: 

Edistanceih j) = E^x-elec x k + Samp X fc X e,;j 2 (8) 

k is the size of the packet. E es ti ma te(h j) denotes the energy of node j 
estimated by node i, which can be estimated according to: 



771 /• -\ 771 ^init E estimate {i t j) r™. /. -\ , n /n\ 

E estimated, J) = E init — -— x [Times{i,]) + 1 (9) 

1 imes(i,j) 

When a node receives a packet, it will evaluate the remaining energy of all 
the neighbors and updates all the values in the routing table to dynamically 
select the next hop. 

When round is multiple of roundToUpdate, procedure of adjusting the 
pheromones starts. First, the pheromones should be evaporated. We utilize 
the equation (|10[) to evaporate pheromones. 



p stands for the fraction of pheromones that are not evaporated. Then 
the procedure of depositing pheromones is performed. Each node selects the 
neighbor with the maximum estimation energy (e.g. the node j), and increases 
the pheromone of node j by E distance (i, j) as: 

= V(hj) + Edistance(i,j) (H) 

When round is multiple of roundToUpdate, algorithm [3] is called to ad- 
just the pheromones for each node. Firstly, the evaporating and depositing of 
pheromones are taken(i.e. lincQ]- line [9]). Then updating of the routing table 
are carried out (i.e. line fTTI - line fT4|) . 

Once algorithm [3] is finished, in the next period of roundToUpdate, the node 
with the highest E distanced j) will have a higher probability of being selected 
as the next hop. After finishing depositing pheromones, the process of updating 
the estimated energy value is performed. Each node will broadcast its current 
energy in the range of R. The value of E est i ma te(i, j) will be updated. 

If a node happens to exist in the conjunction of two or more different paths, 
the parent nodes will also use equation to increase the pheromones. 

Basic-DAACA algorithm avoids the computing overhead of repeatedly con- 
structing or maintaining the network topology, which saves energy and prolongs 
the network lifetime. In the following sections, some heuristic methods are 
proposed to enhance the performance of Basic-DAACA. 
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Algorithm 3 Adjusting pheromones of Basic-DAACA 
l: for all i e V do 

2: i evaporates the pheromones according to equation (|T0|) . 
3: i searches for the node with the highest E es timate{i, j) where j € Nb(i) in 
neighbor set. 

4: i deposits pheromones according to equation l[TT]l . 
5: if conj(i) ^ 2 then 
6: for all t G Nb(i) do 

7: t deposits pheromones of i according to equation 

8: i broadcasts the current energy. 

9: end for 
10: end if 

11: if i receives a broadcast message from j then 

12: £ updates E est i ma te (hj)- 

13: £ updates the routing table. 
14: end if 
15: end for 



5.3. Elitist Strategy based DAACA(ES-DAACA) 

In order to establish better aggregation paths, the global pheromones de- 
position is employed. We call this heuristic method Elitist Strategy based 
DAACA(ES-DAACA for short). We apply another format of the packet header, 
which obtains the id sequence of nodes and total energy consumption as: ID -List, 
Energy-Consumption. ID-List records all the nodes in the trace from the 
source to the current node. Energy-Consumption is defined as the total trans- 
mission energy costed from the source to the current node. 

Once a packet is received by a node, algorithm 0] will be taken which adds id 
information and energy consumption information into the newly defined packet 
header. Owing to the header, the sink node can find the minimum cost ag- 
gregation path and adjust the amount of pheromones to obtain more optimal 
solutions. 

Algorithm 4 Packet receiving process of ES-DAACA 
1: if £ receives the data packet then 
2: i selects the next hop j. 
3: i adds its id information into ID-List. 
4: i adds Edi sta nce{i, j) into the energy consumption field. 
5: end if 



When a packet is sent to the sink node. The sink node checks if the 
Energy -Co sumption is smaller than Mpcost, if so, the Mpcost and the corre- 
sponding ID-List will be updated. When round is multiple of roundToUpdate, 
the sink node will perform global deposition process. Each node in the ID-List 
will add additional Mpcost pheromones as in equation ([12]) , Therefore, ES- 
DAACA can increase the amount of phermones of those nodes who locate in 
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the path with the minimum energy consumption. More optimal solutions of 
data aggregation paths will be likely to be built. Obviously ES-DAACA will 
enhance the performance of Basic-DAACA. 

Vihj) = Vihj) + Mpcostij G Nb(i),i G PminJ £ Pmin) (12) 

Therefore, ES-DAACA will enhance the probability of finding global opti- 
mization solution, which preserves the energy and prolongs the network lifetime. 

5.4. Maximum & Minimum based DA AC A (MM-DAACA) 

To limit the bound of the pheromones, we impose two constants to represent 
the maximum and minimum bounds of pheromones. We name this heuristic 
algorithm as Maximum & Minimum based DAACA (MM-DAACA for short). 

In MM-DAACA, all the steps are identical with ES-DAACA, the only dif- 
ference is that we use [r]Min, VMax] to limit the range of the pheromones. The 
reason is that if the pheromone is not limited in a range, some paths will own 
higher probabilities than the others, nevertheless, with the transmission going 
on, nodes in this kind of path will cost more energy than the others. But accord- 
ing to equation ([3]), the amount of pheromones is still large, which may cause 
them more likely to be selected as the next hop, and ultimately results in local 
optimal solution. However, if the constraints of pheromones are imposed, the 
aforementioned phenomenon can be avoided, variety of paths will be formed, 
hence, the global optimal solution will be more likely to be found. 

5.5. Ant Colony System based DAACA(ACS-DAACA) 

Similar to the ant colony system mentioned in section 13. 2\ we develop the 
Ant Colony System based DAACA (ACS-DAACA for short) algorithm for es- 
tablishing the data aggregation paths to preserve energy. 

In ACS-DAACA, when selecting the next hop node, the algorithm [5] is per- 
formed. q is a parameter. After successfully sending the packet in each round, 



Algorithm 5 Selecting the next hop of ACS-DAACA 
1: Generate a random number q in [0, 1]. 
2: if q ^ go then 

3: Select the neighbor node with the maximum r(i,j) Q x ri(i,j) 
4: else 

5: select the next hop according to (1) 
6: end if 



each node will locally update the pheromones as follows. 

r](i,j) = (1 - p) x r)(i,j) + p x MpcostU e Nb(i),i e P m in,3 S P mm ) (13) 



Since the result of equation (| 13[) is proved to be constrained within a certain 
range 16| , therefore, it is similar to MM-DAACA which can maintain the variety 
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of the solutions. What's more, in ACS-DAACA, we propose local pheromone 
updating strategy as : 

V{iJ) = (1 ~ x r)(i,j) + C x Ar)(i,j) best (14) 

£ stands for the rate of evaporating pheromones. Where 

Arj(i,j) best = min (r)(i,j)) (15) 

jeNbr(i) 

The local pheromones updating aims at reducing the probability of the se- 
lected node, thus the probabilities of unselected nodes will increase. Conse- 
quently, ACS-DAACA will generate more aggregation paths than MM-DAACA 
through which global optimal solution will be more likely to be found. 



6. Simulation Results and Analysis 

In this section, several simulations are performed to evaluate and compare 
the performances of the DAACA family with other data aggregation algo- 
rithms. We use Java to develop a simulation platform to analyze the energy 
cost of each node and the network lifetime. We have testified the accuracy 
of the platform by simulating the third set of experiments in [14| and get 
identical statistics on the platform. Then we simulate different algorithms(i.e. 
LMST, PEDAP, PEDAP-PA, L-PEDAP, ACA, Basic-DAACA, ES-DAACA, 
MM-DAACA, ACS-DAACA) on the platform. In the experiments, multiple 
sources periodically send data to the sink node. Related parameters are shown 
in table [3] which are identical with those in [14J . We have run millions of ex- 
periments to highlight the advantages of DAACA in this paper. What's more, 
features of network topology are discussed. Furthermore, the time complexity 
analysis and characteristic analysis are given. 



6.1. Average Energy Cost 

All the nodes are deployed randomly as (Length x Width, Node Number). 
Different numbers of packets are transmitted from the source nodes to the sink 
node. We calculate the average remaining energy of nodes with different algo- 
rithms. Higher remaining energy means less energy cost. Figure [2] - figure [4] 
illustrate the comparison results. Firstly each algorithm will initialize the topol- 
ogy of the network. With the transmissions going on, some of the algorithms 
will re-construct (or maintain) the topology to balance the energy cost of each 



node [17] . Hence, the topologies of some algorithms are dynamical. 



The mechanism of establishing network topology of LMST and L-PEDAP 
are nearly identical, the energy consumptions of construction and maintenance 
of the network topology are high. Because each node in the network establishes 
the minimum spanning tree in the neighborhood locally. Then local minimum 
spanning trees connect into the topology of the network. When establishing 
the local minimum spanning tree, the weight of edge is directly related to the 
remaining energy of the node. Hence, each node needs to communicate with 
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Table 3: Experimental Parameters 



Item 


Parameters 


Algorithm 


LMST, PEDAP, PEDAP-PA, L- 




PEDAP, ACA, Basic-DAACA, ES- 




DAACA, MM-DAACA, ACS-DAACA 


Packet Size 


4098 bits 


(Length x Width, Node Number) 


(40 x 50, 200), (60 x 60, 400), (70 x 




80, 600), (80 x 90, 800), (100 x 100, 




1000), (100 x 120, 1200), (100 x 140, 




1400), (120 x 130, 1600), (130 x 120, 




1800), (140 x 140, 2000) 


Node Distribution 


Random 


Initial Node Energy 


10J 


Packets Number 


1000, 2000, 3000, 4000, 5000, 6000, 




7000, 8000, 9000, 100000 


Number of Sources 


10 


R 


10 


a,P 


a = 2,P = 2 


Initial Pheromone 


0.8 


T]Mim VMin 


TjMin = 0.5, TjMax = 0.9 


P 


p = 0.2 


roundToU pdate 


100, 200 


c ' 


0.9 




0.5 


rho in ACA 


0.3 




Figure 2: Average Energy Cost after Transmitting 10000 Packets 
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Figure 3: Average Energy Cost of the Network with the size (140 X 140, 2000) 




Figure 4: Average Results of all the Experiments of Average Energy Cost 
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neighbors to obtain their remaining energies periodically to update the weights 
of edges locally, frequent communications cast a heavy burden of all the nodes 
in the network. Therefore, the LMST and L-PEDAP cost more energy than 
others. 

In ACA, the energy cost concentrates on the frequent depositions of pheromones. 
Because whenever a packet is received by a node, if the TTL of the packet is 
bigger than the i?iJC(Extended Hop Count) Q of that node, the deposition 
of pheromones will start. The feedback message will be transmitted to inform 
the upstream nodes to deposit pheromones. Suppose there are a large number 
of packets transmitting in the network, and TTL >> EHC, the feedback mes- 
sages will be massively used, which consumes considerable energy. Therefore the 
energy consumption mainly concentrates on too frequent transmitting feedback 
messages. 

In PEDAP and PEDAP-PA, the Prim algorithm are used to construct the 
topology and maintain the topology from the global view. Such methods can 
save the energy of communicating with neighbor nodes, but the sink node will 
suffer from more energy consumptions than the other nodes, because this is a 
centralized version of method. The sink node is responsible for establishing the 
topology of the network and informing all the nodes about the topology of the 
network via broadcasting. The energy cost of broadcast is high, especially when 
some nodes locate far away from the sink node. Whenever the network topology 
is changed, the broadcasting will be executed. Therefore, the energy costs of 
PEDAP and PEDAP-PA in constructing and maintaining the network topology 
are high. 

In the family of DAACA, there is no energy consumption in maintaining the 
topology or reconstructing the MST, therefore, the average energy cost is lower 
than other algorithms. The energy consumption concentrates on sending packets 
to the next hop. After certain rounds of transmissions, the global pheromones 
adjustment will be carried out, which does not cost much energy. From the 
figure [5]- figure 0] we can see, the relationship among the four methods of the 
DAACA family in terms of average energy cost is ACS<MM<ES<Basic. MM- 
DAACA, ES-DAACA and ACS-DAACA enhance the quality of the solutions of 
Basic-DAACA which generate more energy preservation data aggregation paths 
for the network and eventually save the energy. 

6.2. Energy Difference 

In this section, we calculate energy difference to evaluate the balanced con- 
sumption of energy. We aims at testifying which method can balance the usage 
of energy. Because the balanced consumption can directly affect the lifetime of 
the network, which is the most prominent feature in energy efficiency in wireless 
sensor networks. If the energy difference is big, the energy is not averagely used 
whereas smaller energy difference indicates the balanced energy consumption. 

In figure [5] - figure [7J we can see that the energy difference of L-PEDAP is 
the highest, the reason is that in L-PEDAP, each node needs to establish local 
minimum spanning tree periodically. L-PEDAP is energy aware, the weights of 
edges in the neighborhood always change, which demands each node to set up 
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Figure 5: Energy Difference after Transmitting 10000 Packets 




Figure 6: Energy Difference of the Network with the size (140 X 140, 2000) 
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Figure 7: Average Results of all the Experiments of Average Energy Cost Energy Difference 



the local minimum spanning tree repeatedly. Although L-PEDAP can balance 
the consumption of energy within neighborhood, a node with bigger neighbor set 
will actually consumes more energy than the smaller one. Because in the dense 
region of the network, each node owns many neighbor nodes. When establishing 
the LMST, each node needs to communicate with all the neighbors to acquire 
the energy information of the neighbors. Therefore, a node with large neighbor 
node set will spend much energy in establishing LMST, whereas a node owns 
small neighbor set cannot spend much energy, which yields unbalanced energy 
consumption such two kinds of nodes. Therefore, a big energy difference occurs. 

Although LMST utilize RNG to establish the network topology, the same 
problem of L-PEDAP still exists, which leads to a big energy difference. 

In PEDAP and PEDAP-PA, the sink node consumes most of the energy, 
because it needs to broadcast the topology of the network to all nodes in the 
network periodically. Thereby, there is a large energy difference between the 
sink node and others. 

Whereas in DAACA, it takes the advantages of both centralized and decen- 
tralized algorithms. When selecting the next hop, each node will estimate the 
remaining energy of its neighbor to achieve balanced usage of energy. Moreover, 
in the process of depositing global pheromones, the path with minimum energy 
cost will obtain additional pheromones, by which the balanced usage of energy 
is implemented globally. The results show that DAACA can balance the energy 
cost of the network. The energy difference relationship among DAACA family 
is Basic>MM>ES>ACS. 
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LMST PEDAP PEDAP_PA L_PEDAP ACA Basic_DAACA ES_DAACA MM_DAACA ACS_DAACA 

Figure 8: Lifetime of Network 

6.3. Network Lifetime 

The initial energy of each node is set as 0.05 J. We measure how many rounds 
the network will sustain until any node exhausts in each algorithm. The pa- 
rameter roundToUpdate is set as 100 and 200 respectively. From figure [U the 
DAACA family shows longer lifetime than the other algorithms. And the rela- 
tionship within the DAACA family is ACS>MM>ES>Basic. In PEDAP and 
PEDAP-PA, the sink node will firstly exhausted because of frequent broadcast- 
ing of topology. The lifetime of L-PEDAP and LMST are short which mainly 
results from massive energy costs in establishing local minimum spanning tree. 
In ACA, the energy cost mainly concentrates on some shortest paths. DAACA 
can guarantee the balanced energy costs of nodes by adjusting quantities of 
pheromones, therefore, it can prolong the lifetime of the network. 



6.4- Successful Ratio of One Hop Transmission 

In this section, we come to analyze the average success ratio of one hop 
transmission of each algorithm. We choose the data set (100 x 100, 1000). 
Packets are sent from the sources to the destination nodes. We count the success 
ratio of one hop transmission of each packet and calculate the average successful 
ratio of the one hop transmission. In figure [9l the successful ratio of DAACA 
is higher than other algorithms. The relationship within the DAACA family is 
ACS>MM>ES>Basic. 



6.5. Network Topology Analysis 

In this section, we analyze the performance of network topology in terms 
of node average degree and average transmission radius. As the experiment in 



12j, the data set (100 x 100, 1000) is chosen to run the simulation. 



6.5.1. Average Degree 

In this experiment, we come to analyze the average degree of each node. The 
average degree is calculated as the number of edges divided by the number of 
nodes. Bigger average degree indicates the higher spatial resuepjj- 

Table [4] lists the average degree of each node in the network. In PEDAP and 
PEDAP-PA, the Prim algorithm is used to generate the minimum spanning 
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Figure 9: Average Successful Ratio of One Hop Transmission 



Table 4: Average Degree 

Algorithm Average Degree 

LMST 2.04 

PEDAP & PEDAP-PA 0.999 

L-PEDAP 2.04 

ACA 0.999 

DAACA 0.999 
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tree of the network globally, since in the minimum spanning tree, the number 
of edge equals to the number of nodes minus 1, therefore the sum of the degree 
equals to the number of nodes minus 1 and the average degree is n — 1/n where 
n is the number of nodes in the network. In ACA and DAACA, although the 
topologies are dynamic, each time only one node is selected as the next hop until 
the packet arrives at the sink node. Therefore the total degree still equals to 
that of PEDAP and PEDAP-PA. Whereas in LMST and L-PEDAP, local MST 
is produced, which cause many loops in the topology, in such circumstance the 
number of edges is bigger than n. Hence, the average degree is higher than 
other algorithms. 

In a word DAACA can achieve a higher spatial reuse. 



6.5.2. Average Transmission Radius 

In this experiment, we calculate the average transmission radius of each 
node. As illustrated in 12[, average transmission radius is another factor which 
indicates the spatial reuse. 



Table 5: Average Transmission Radius 



Algorithm Radius 

LMST 4.81 
PEDAP & PEDAP-PA 2.19 
L-PEDAP 2.07 
ACA 1.8 
DAACA 1.79 



As illustrated in 12| , the smaller average transmission radius will lead to a 



higher spatial reuse. In table El the results of the average transmission radius 
of each algorithm are listed. In LMST and L-PEDAP, the local MST are es- 
tablished, more links will be used to establish the topology, which is the main 
reason of causing the higher transmission radius. In PEDAP and L-PEDAP, the 
average degree is 2.19 and 2.07 respectively. In ACA and DAACA, the next hop 
node is selected which is nearer to the sink node than to the current node. The 
selection is based on the probabilities; the shorter edge has a higher possibility 
to be chosen. In DAACA the edge is selected based on distance and energy, 
therefore, the radius is a little shorter than ACA. 

From the topology analysis, we can conclude that DAACA achieves a higher 
spatial reuse than the other algorithms. 



6.6. Time Complexity Analysis 

In this section, we mainly analyze the complexity of constructing and main- 
taining the topology. Because in the initialization phase, the steps (e.g. initialize 
the neighbor set, initialize routing protocol, et al.) are nearly identical, we only 
compare the time complexity in constructing and maintaining the topology of 
each method. 
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6.6.1. Time Complexity of Topology Construction 



Table 6: Complexity of Topology Construction 

Algorithm Time Complexity 

LMST 0(nelog(n)) 

PEDAP-PA k LEDAP-PA 0{elog(n)) 

L-PEDAP 0{nelog{n)) 

AC A 0{n 2 ) 

DAACA o\n 2 ) 



Table [6] reveals the time complexity of each algorithm in constructing net- 
work topology. In LMST, each node establishes the local minimum spanning 
tree, the Prim algorithm is utilized to generate the MST, the complexity to 
set up MST is 0(elog(n)), where e is the number of links within the visible 
neighborhood [l3[ and e > n. But LMST is only apt to the network with a 
small number of nodes or low density network. Because when LMST is taken 
upon numerous nodes, each node will own a large set of neighbors which yields 
to e ps n 2 . Therefore the overall complexity in constructing LMST is high. 
Whereas in PEDAP and PEDAP-PA, the centralized version of minimum span- 
ning tree is established, which owns a complexity of 0(elog(n)). Compared with 
LMST, the complexity is relatively low, but it is a centralized method which 
casts a heavy burden on the sink node. In L-PEDAP, the topology construc- 
tion is identical with LMST; therefore, the complexity is still high. In ACA, 
when establishing the network topology, each node must broadcast messages 
to his neighbors in order to inform them its hops to the sink node. Therefore, 
the complexity is 0{n 2 ). Topology construction DAACA is nearly the same as 
ACA, thus the complexity is 0(n 2 ). 



6.6.2. Time Complexity of Topology Maintenance 

In this section, the complexity of maintaining the topology is discussed. Ta- 
ble [7] lists the complexity of each algorithm. Since each method updates its 
topology in every roundeToUpdate rounds periodically, it is feasible to com- 
pare the complexity in just one period. In LMST, PEDAP, PEDAP-PA and 
L-PEDAP, all the topologies are required to be re-constructed. Thus the com- 
plexities are the identical with those in constructing topologies of section [6. 6. II 
In ACA, whenever a node meets the condition of EHC < TTL, operations 
of pheromones are performed, some extra control messages [14j are needed 
to be transferred, therefore, the complexity of ACA is 0(rn) (r stands for 
roundTolIpdate). Whereas in the family of DAACA, each node locally cal- 
culates the probability of selecting the next hop in every round. Only one 
broadcasting message for updating the energy evaluation is required, therefore, 
the complexity is 0(n). From table [3 it can be concluded that the complexity 
of maintaining the topology of DAACA is lower than the other methods. 
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6.7. Other Characteristic 

In this section, we analyze other characteristics of DAACA. 

6.7.1. Robustness 

When selecting the next hop, each node will refer to the probability recorded 
in the routing table, all the path are dynamical, if a node is removed from 
the network, only the nodes nearby are required to update the routing table. 
Therefore, the adjustment of removing a node can be achieved locally, therefore, 
DAACA owns the feature of robustness. 

6.7.2. Fault Tolerance 

Suppose errors occur in the element of the node which prohibits it from 
working. In this case, the node nearby can remove this node from the routing 
table. Therefore, the faulty nodes will no longer exist in the topology of the 
network. Hence they cannot affect the transmission of the network. And we 
can conclude that DAACA is fault tolerant to faulty nodes. 

6. 7. 3. Scalability 

Since nodes are deployed randomly in the network, and each node only needs 
to maintain the routing table for dynamically selecting the next hop. Therefore, 
if a node is added in the network, it only needs to broadcast his identity infor- 
mation and sets up routing table, the nodes nearby are only required to add an 
new item of the routing table. So we can conclude that DAACA is scalable. 

7. Conclusion 

In this paper, a family of ant colony algorithms called DAACA for data 
aggregation has been presented which contains three phases: the initialization, 
packet transmission and operations on pheromones. After initialization, each 
node estimates the remaining energy and the amount of pheromones to com- 
pute the probabilities used for dynamically selecting the next hop. After certain 
rounds of transmissions, the pheromones adjustment is performed periodically, 
which combines the advantages of both global and local pheromones adjustment 
for evaporating or depositing pheromones. Four different pheromones adjust- 
ment strategies are designed to achieve the global optimal network lifetime, 



Table 7: Complexity of Topology Maintenance 

Algorithm Time Complexity 

LMST 0{nelog(n)) 

PEDAP-PA & LEDAP-PA 0(elog(n)) 

L-PEDAP 0{nelog{n)) 

AC A 0(rn) 

DAACA 0{n) 
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namely Basic-DAACA, ES-DAACA, MM-DAACA and ACS-DAACA. Com- 
pared with some other data aggregation algorithms, DAACA shows higher su- 
periority on average degree of nodes, energy efficiency, prolonging the network 
lifetime, computation complexity and success ratio of one hop transmission. At 
last we analyze the characteristic of DAACA in the aspects of robustness, fault 
tolerance and scalability. 
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